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ABSTRACT 

Aims. We present new millimetre 43 GHz observations of a sample of radio-bright Planetary Nebulae. Such observations were carried out to 
have a good determination of the high-frequency radio spectra of the sample in order to evaluate, together with far-IR measurements (IRAS), 
the fluxes emitted by the selected source in the millimetre and sub-millimetre band. This spectral range, even very important to constraint the 
physics of circumstellar environment, is still far to be completely exploited. 

Methods. To estimate the millimetre and sub-millimetre fluxes, we extrapolated and summed together the ionized gas (free-free radio emission) 
and dust (thermal emission) contributions at this frequency range. By comparison of the derived flux densities to the foreseen sensitivity we 
investigate the possible detection of such source for all the channels of the forthcoming ESA's PLANCK mission. 

Results. We conclude that almost 80% of our sample will be detected by PLANCK, with the higher detection rate in the higher frequency 
channels, where there is a good combination of brighter intrinsic flux from the sources and reduced extended Galactic foregrounds contam- 
ination despite a worst instrumental sensitivity. From the new 43 GHz, combined with single-dish 5 GHz observations from the literature, 
we derive radio spectral indexes, which are consistent with optically thin free-free nebula. This result indicates that the high frequency radio 
spectrum of our sample sources is dominated by thermal free-free and other emission, if present, are negligible. 



- Planetary Nebulae: general radio continuum: stars 



Key words, stars: Planetary Nebulae 

1. Introduction 

The PLANCK ESA mission will provide us with nearly full 
sky maps over a wide range of frequencies, from 30 to 900 
GHz. Therefore, even if designed for cosmological studies, the 
mission will have profund impact on fundamental Physics and 
Galactic and extragalactic astrophysics. Planck will be sensi- 
tive to the millimetre emission from dusty envelopes of stars 
and we expect the dusty circumstellar envelopes, that char- 
acterize the latest stages of stellar evolution, to be source of 
relevant foreground contamination. When a low or intermedi- 
ate mass star is approaching the end of its evolution, it goes 
through a period of heavy mass-loss known as Asymptotic 
Giant Branch (AGB) phase. The ejected envelopes are par- 
tially condensed in dust grains and completely obscure the cen- 
tral star. Immediately after the AGB evolutionary phase, the 
mass-loss stops and the stars may become optically visible as 
the dusty shells disperse (Proto Planetary Nebula, PPN phase). 
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Eventually, once it reaches a temperature of 2 - 3 xlO K, 
the central star starts to ionize the AGB shell and a Planetary 
Nebula will form. Dusty envelopes re-radiate the absorbed stel- 
lar light showing a clear signature in the far-infrared spectrum, 
i.e. an IR excess with peculiar IRAS colours. In addition to that, 
PNe show a radio continuum due to free-free emission from the 
fraction of the CSE ionized by the central star. 

From a feasibility study, Umana et al. (2006) concluded that 
a sizable 300) sample of AGB and post-AGB stars would 
be detected during the mission and derived estimates for the 
expected flux densities at various Planck channels. However, 
the simulations carried out by Umana et al. ( I2006l l, on a sub- 
sample of PNe rely only on NVSS fluxes, obtained at 1.4 GHz, 
extrapolated to the Planck frequencies. This leads to underesti- 
mate the contribution due to free-free as, at this frequency, PNe 
are often optically thick (Siodmiank & Tylenda 2001, Luo et 
al. l20Q5l 

PNe are among the brightest Galactic radio sources. Some 
of them could also reach a flux density of Jy level. More than 
800 have been detected at least at one frequency and for 200 
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Table 1. The selected sample 
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morphological and spectral information (between 1.4 and 22 
GHz) were obtained. Most of the higher frequency (22 GHz) 
data were obtained with interferometers (i.e. VLA) while very 
little single-dish, high frequency measurements are available. 

In this paper we present new 43 GHz, single-dish observa- 
tions by using the 32 m INAF-IRA Radiotelescope at Noto of a 
sample of PNe, which are potential foregrounds for PLANCK. 
The main goal of this project is to obtain reliable estimates of 
flux density expected at PLANCK channels, by building and 
modeling their spectral energy distribution (SED). This in turn 
will contribute to the compilation of the PLANCK pre-launch 
catalogue. We stress here that 43 GHz is one of the observing 
channels of the forthcoming PLANCK mission. Therefore, at 
this frequency band, we would obtain a direct measurement of 
the expected flux and not an extrapolation. 

As added value, the present work provides the first size- 
able dataset of 43 GHz measurements of PNe, that constitute 
strong constraints to the observed SEDs in the very impor- 
tant spectral region where free-free emission and thermal dust 
emission may overlap. While interferometric high frequency 
observations provide us with detailed morphological informa- 



tion they quite often fail to entirely recover the extended emis- 
sion. This eventually leads to underestimate the total radio flux 
density. This problem is overcame by single-dish observations. 
Assesting the SEDs in the radio-millimetre spectral range is 
a crucial step for the study of the physics of dusty envelopes 
around PN. A correct evaluation of the free-free contribution, 
up to millimetre range, when combined with information pro- 
vide by far-IR observations, would allow to determine the pres- 
ence of an excess due to the presence of a cold dust compo- 
nent/s (Gee et al. [19841 Hoare et al. [19921 Kemper et al. |20"0"2| | 
or of alternative emission mechanisms (Casassus et al.[2007). 
Since PNe and their progenitors are believed to be among the 
major sources of recycled interstellar matter, determining the 
properties of the dust ejected in the ISM is very important to 
the study of the Galaxy evolution in general. 

2. The 43 Ghz Noto Survey 

2.1. Sample Selection 

Our sample has been selected mainly from Condon & 
Kaplan ( 1998), who performed a cross-correlation between the 
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Table 2. Results 
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* References for 5 GHz measurements: A) Gregory et al. (19961; B) Griffith et al. fl995l ; C) Griffith et al. ( fT994| l; D) Wright et al. ([1994); 
E) Becker et al. ( fl99"lT >; F) Haynes et al. ( [19791 ; G) Milne ( 1979 ); H) Milne ( fl975T > 



Strasbourg-ESO Catalogue of Galactic Planetary Nebulae (Cat. 
<V/84>, Acker et al. [1992j l and the 1.4 GHz NRAO VLA Sky 
Survey (NVSS). Among these sources we selected only those 
whose flux densities at 1.4 GHz (NVSS) is higher than 100 
mJy, for a total of 64 PNe. In the conservative hypothesis of 
optically thin nebula at 1.4 GHz, a cut-off at 100 mJy guar- 
antees a 43 GHz flux density easily detectable with the Noto 
Radiotelescope. This estimated 43 GHz flux density will be a 
lower limit in the case of optically thick nebula at 1.4 GHz. 
We also considered Condon et al. ( 1999 ) who selected infrared 
PNs in NVSS by cross-correlation between NVSS and a sam- 
ple from IRAS PSC on the basis of infrared colours characteris- 
tic of PNe. Only 3 over the 122 of infrared PNs candidate have 
S 1.4 GHz ^ 100 mJy, yielding to a final sample of 67 objects. 

To avoid problems due to possible contamination in the 
Noto beam (HPBW= 54"), for each source of the sample, we 
have extract from NVSS a 25' x 25' field, centred at the po- 
sition of the target. Sources that show a very extended (> 2') 
emission or that are located in very high confusion region have 
been rejected. This reduces our sample to 62 objects. The list 



of the selected targets, with names and positions, is reported in 
TableQ] 

2.2. Observations and Results 

In the last few years the INAF-IRA 32m Noto radiotelescope 
has been subjected to a series of structural improvements which 
have increased remarkably its potential capabilities as single 
dish instrument. In particular, the installation of an active sur- 
face allows to operate with good performances at high frequen- 
cies. 

The observations reported in this paper were carried out 
in different epochs, between 2005 and 2006. The 43 GHz su- 
pereterodyne receiver is cooled double polarization receiver, 
with typical zenith system temperature (T sys ), in both channels, 
of the order of 80 - 100 K, depending on the weather condi- 
tions. The gain ranges from 0.05 to 0.1 K/Jy (Leto et al. 2006), 
depending on elevation, and this determines a zenith System 
Equivalent Flux Density (SEFD) of 1600 - 2000 Jy. The obser- 
vations were performed with a 400 MHz instantaneous band. 
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All the sources of our sample were observed with the on the 
fly (OTF) scan technique, which consists in driving the beam 
of the telescope across the source in RA direction. The typical 
scan duration was of the order of 20 sees, short enough to re- 
main close enough to the white noise regime of the radiometer. 
In order to achieve a good signal to noise ratio, each source was 
observed many times, for a total integration time of 30 minutes. 
Multiple OTF scans were then added together reaching a typi- 
cal rms of 2 - 3 mK. 

Daily gain curves were obtained and the flux scale was 
fixed by using NGC7027 as primary calibrator and 3C286 as 
secondary calibrator. The adopted flux densities for NGC7027 
and 3C286 are respectively: 5.07 and 1.86 Jy. The assumed 
flux densities are those reported by Ott et al. (1994); in the 
case of NGC7027 the measure has been corrected for the ob- 
served decreasing of 0.15 percent/year as reported by Perley et 
al. (120061 . 

We detected 55 out of the 62 observed objects, with a 89 
% detection rate. Results of such observations are reported in 
Table [2j where the 43 GHz measured flux density, or its 3<x 
upper limit, is listed. In the 4 th and 10 th columns of the same 
table we report the angular size at 1.4 GHz of each source 
(01.4 GHz), as derived from the analysis of the relative NVSS 
map, which has an angular resolution of 45". The source angu- 
lar size has been derived as the geometrical mean of the minor 
and major axes obtained by fitting a two dimensional Gaussian 
at the source position in the map of the brightness distribution, 
by using the task JMFIT of the NRAO Astronomical Image 
Processing System (AIPS). 

To take into account for possible partial resolution of the 
source by the Noto beam, assuming that the source size at 43 
GHz is the same as measured at 1.4 GHz, we corrected the 
measured flux density as 



S c 43 GHz - S 43 GHz X 



^Noto + ^1.4 GHz 
^Noto 



where #1.4 ghz and #Noto are the width of the source at 1.4 
GHz (reported in Table |2| and the width of the Noto radiote- 
lescope beam (HPBW=54") respectively. Sources having an- 
gular size greather than 25", comparable to the 43 GHz Noto 
beam (HPBW), have flux corrections significantly larger than 
the rms associated to the flux density. This means that they can- 
not be regarded as point-like. The resulting fluxes S c 43 ghz are 
listed in Table [2] 
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Fig. 1. Calculated spectral index as function of flux at 43 GHz. 
As expected major dispersion is evident at lower fluxes 



available only for 41 targets. The 5 GHz single dish measure- 
ments are from surveys performed by NRAO Green Bank 91 
cm telescope, with a beamwidth (HPBW) * 3.5' (Gregory et 
al. [1996 Becker et al. 1991) and Parkers 64 m telescope with 
a beamwidth « 4.5' (Griffith et al. [19951 Griffith et al. [T9941 
Wright et al. [T994l Haynes et al. [19791 Milne [19791 and Milme 
& Aller l 19751 ). as indicated in Table[2] For the sources observed 
in more than one survey, we used the average of the fluxes from 
the different measurements. 

The spectral indices a, calculated between the two frequen- 
cies, are shown in Fig.Q]as a function of the flux density at 43 
GHz. Since its accuracy depends on the value of the measured 
flux densities and associated uncertainties, we computed the 
expected 1 <x as a function of the flux density, assuming a typi- 
cal rms, associated to the 43 GHz and 5 GHz measurements, of 
o"43GHz ~ 50 mjy and a <t 5 ghz ~ 35 mjy, respectively. The +1 cr 
uncertainty values around the mean value a = -0.09 are shown 
in Fig. [T]as dot-dashed lines. For about 70% of the sources we 
get a spectral index inside the uncertainty lines, showing that it 
is statistically consistent with the value a = -0.1, typical for 
an optically thin free-free emission. We may conclude that for 
our sample, at least up to 43 GHz, the SED is dominated by 
the free-free emission and other contributions, if present, are 
negligible. 



3. The SEDs 

3.1. The free-free contribution 

As a first step in our analysis, we derived the spectral index 
a, (S v oc v"), by combining the fluxes measured at 43 GHz 
with the 5 GHz single dish measurements available from liter- 
ature. In order to prevent any error due to resolving out some 
of the flux density, we considered only the objects that have an- 
gular size, at 1.4 GHz (G\a ghzX lower than 25". This reduces 
our sample detected sources to 42 PNe that can be, quite con- 
fidentially, considered point-like with respect to the telescope 
beam at the observational frequency. Literature 5 GHz data are 



3.2. The dust contribution 

PNs are usually surrounded by a dusty envelope which is the 
remnant of the precursor's wind. Dust thermally re-radiates the 
UV radiation of the central star determining an excess that may 
extend from the far-IR to the radio region. Such a contribution 
is typically of the order of 40% of the total flux from a PNe 
(Zhang & Kwok 1991 ), being more important in young PNe, 
since, in more evolved PNe, the circumstellar material has al- 
ready dispersed. We built the dust component of our sample of 
PNe by using IRAS measurements. 

The IRAS data have been fitted by assuming that the dust 
emits as a blackbody modified by the frequency dependent dust 
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Fig. 2. Spectral Energy Distributions of four of our targets. 
Both ionized gas (dotted line) and thermal dust (dashed line) 
contributions are shown. Continuous line is the total emission. 
Diamonds: IRAS data; squares: our measures at 43 GHz (7 x 
10 3 pm); triangle: literature measures at 5 GHz (6 x 10 4 pm); 
cross: NVSS data at 1.4 GHz (2 x 10 5 pm) 

opacity, that is F v oc v p B v (T$), where B v {Td) is the Planck func- 
tion for dust temperature 7d, and p is the emissivity index. The 
index p strongly depends on the mineralogical composition of 
the grain and on their physical shape. In our calculation, fol- 
lowing the results from detailed modelling of dust emission 
from PNe (Hoare 1990 ), we assumed, for the slope of the far-IR 
emissivity law, a typical value of p- 1.5. Because strong con- 
tamination by line emission of 12 pm IRAS band (Stasinska & 
Szczerba 1999) only 25, 60 and 100 pm fluxes have been used 
in the fitting procedure. 

4. Estimates of flux density at Planck channels 

The SEDs with both ionized gas and dust contribution, evalu- 
ated in Sect. 13.11 and 13.21 were built for all the targets in our 
sample. In Fig. [2] four of such SEDs are shown as an exam- 
ple. In the SEDs, between radio and far-IR there is a large 
gap where measurements are missing. New millimetre and sub- 
millimetre observations are clearly necessary to better charac- 
terize and constraint the emission from different components, 
where the radio flux, from the ionized fraction of the nebula, 
could be merged to the dust contribution. This gap will be par- 
tially filled by the ESA-PLANCK mission, as the satellite is 
equipped with a Low Frequency Instruments (LFI), operating 
at 30, 44 and 70 GHz, and a High Frequency Instruments (HFI), 
operating at 100, 143, 217, 353, 545 and 857 GHz. 

In order to evaluate the possibility that PLANCK will ac- 
tually detect sources in our sample, we should compare, for 
each observational band, the expected flux from the sources 
to the foreseen sensitivity. The expected fluxes are derived by 
summing, at each PLANCK channel, both dust and ionized gas 
contribution. The first has been obtained by extrapolating the 
modified blackbody fit (dashed lines in Fig. the latter has 
been obtained from fitting free-free emission from an optically 
thin nebula to the 5 GHz (6 x 10 4 pm) and 43 GHz (7 x 10 3 pm) 



Fig. 3. Histogram showing the fraction of planetary nebulae be- 
longing to our sample that can be detected by the Planck satel- 
lite 

data points and then extrapolating to the Planck channels (dot- 
ted line in Fig.EJ, the NVSS measures at 1.4 GHz (2 x 10 5 pm) 
has been also displayed. 

Following Umana et al. (2006), the foreseen total sensitiv- 
ity is assumed to be the sum in quadrature of all the sources of 
confusion noise, namely: the nominal instrumental Planck sen- 
sitivity per resolution element; the Galactic and extragalactic 
foregrounds confusion noise (Toffolatti et al.[l998); the CMB 
confusion noise (Bennett et al. 2003). All the considered con- 
fusion sources will contribute to a total sensitivity (cr) rang- 
ing from 270 mJy for the 30 GHz channel to 220 mJy for the 
857 GHz channel. In Fig. [3] we show the percentage of the 
objects from our sample that could be detected over the 3cr 
threshold for each channels of the PLANCK instruments. A 
small number of PNe belonging to our sample will be detected 
by PLANCK in the radio (LFI channels), while most of our 
targets will be easily seen by PLANCK in the millimetre and 
sub-millimetre bands (HFI channels). We may conclude that 
we foresee an important contribution of the PLANCK mission 
also to the study of PNe, as it would provide measurements in a 
very important, but still far to be fully explored, spectral region. 

5. Summary 

We have presented new 7 mm (43 GHz) observations of a sam- 
ple of radio-bright PNe, carried out with the INAF-IRA Noto 
Radiotelescope. Such observations have been used to derive the 
high-frequency free-free contribution, due to the ionized frac- 
tion of the circumstellar envelope. So far, the majority of high 
frequency radio observations of PNe have been conducted with 
interferometers, that reveal details of source radio morphol- 
ogy but could, in principle, resolve out some of the extended 
emission. Our single dish measurements provide an extended 
database of millimetre observations of PNe, to be used when 
is necessary to know the overall emission, such as when build- 
ing a SED. We used our measurements, that trace the free-free 
contribution, together with IRAS measurements, which trace 
the thermal dust emission, to build up the observed SED, from 
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radio to far-IR and by extrapolation, to estimate the expected 
fluxes in the spectral range between radio and sub-mm, where 
observational data are missing. 

When comparing the expected millimetre-sub-millimetre 
fluxes with the total foreseen sensitivity of the forthcoming 
ESA mission PLANCK we estimate that a consistent number 
of our targets will be easily detected by PLANCK, mostly at 
higher frequency channels. Therefore, even if designed for cos- 
mological study, PLANCK could also contribute to the PNe 
science. Results from such kinds of observations, once mod- 
elled with appropriate code (i.e. CLOUDY, DUSTY), would 
provide important constraints on the chemical composition 
and structure of the CSEs. It would point out the presence of 
multi-shells, related to multiple mass-loss event suffered from 
the central object during its previous evolution (AGB), or a 
contribution of different emission processes, besides free-free 
and thermal from dust, as recently claimed by Cassassus et 
al. ([2007). Moreover, PLANCK results can be considered as 
pathfinder for other future instrumentations, with the same fre- 
quency coverage, such as ALMA, as they will help in planning 
more focused experiments aimed to investigate the morpholog- 
ical details of the sources. 
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